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ABSTRACT 

Aims. The planetary system around the M star Gliese 581 consists of a hot Neptune (Gl 581b) and two super-Earths (Gl 581c and Gl 
58 Id). The habitability of this system with respect to the super-Earths is investigated following a concept that studies the long-term 
possibility of photosynthetic biomass production on a dynamically active planet. 

Methods. A thermal evolution model for a super-Earth is used to calculate the sources and sinks of atmospheric carbon dioxide. The 
habitable zone is determined by the limits of photosynthetic life on the planetary surface. Models with different ratios of land / ocean 
coverage are investigated. 

Results. The super-Earth Gl 581c is clearly outside the habitable zone, since it is too close to the star In contrast, Gl 58 Id is a tidally 
locked habitable super-Earth near the outer edge of the habitable zone. Despite the adverse conditions on this planet, at least some 
primitive forms of life may be able to exist on its surface. Therefore, Gl 581d is an interesting target for the planned TPF/Darwin 
missions to search for biomarkers in planetary atmospheres. 

Key words, stars; individual; Gl 581 — stars; planetary systems — astrobiology 



1. Introduction 

Planets have now been observed around more than 200 main- 
sequence stars. Based on the available observational tech- 
niques, most detected objects are giant (Jupiter-like) planets. 
Therefore, until very recently the existence and possible hab- 
itability of Earth-like planets in extrasolar planetary systems 
was hi ghly speculativ e (Cuntzet al. 2003; Franck et al. 20031; 
von Bloh et al. 2003at iJones et al.] 120(361; [von Bloh et al.l 120071) . 
Bonfils et ah (i2005[) reported the detection of a Neptune size 
planet around Gl 581, an M dwarf star at a distance of 6.26 
pc with a mass of . 3 IMp and a luminosity of 0.013 L©. Very 
recently, lUdrv et all ( |2007|) announced the detection of two so 
called "super-Earth" planets in this system, Gl 581c with a mass 
of 5.06 Me with a semi-major axis of 0.073 AU, and Gl 58 Id 
with 8.3 Me and 0.25 AU. Both mass estimates are minimum 
masses uncorrected for the inclination term sin ;, which is cur- 
rently unknown. 

According to [Valencia et al.l (l2006l) . super-Earths are rocky 
planets from one to ten Earth masses with the same chemi- 
cal and mineral composition as the Earth. In the following, we 
adopt the hypothesis that this is indeed the case, and consider for 
these planets our model previously developed for the Earth, us- 
ing when appropriate scaling laws. This justifies the term super- 
Earth used here for these planets. 

The main question is whether any of the two super-Earths 
around Gl 581 can harbour life, i.e., that any of the planets 
lies within the ha bitable zone (HZ). As a first approximation, 
lUdrv et akl (l2007h computed an equilibrium surface temperature 
for Gl 581c of 20°C for an albedo of 0.5. They neglected, how- 
ever, the likely greenhouse effect of the atmosphere. Typically, 
stellar HZs are defined as regions around the central star, where 
the physical conditions are favourable for liquid water to exist at 



Send offprint requests to: W. von Bloh 



the planet's surface for a period of tiine long enough for biolog- 
ical evolution to occur. iKasting et alj (1 19931) calculated the HZ 
boundaries for the luminosity and effective temperature of the 
present Sun as Ri^ = 0.82 AU and Rqm - 1-62 AU. They defined 
the HZ of an Earth-like planet as the region where Uquid water 
is present at the surface. 

According to this definition, the inner boundary of the HZ 
is determined by the loss of water via photolysis and hydrogen 
escape. The outer boundary of the HZ is determined by the con- 
densation of CO2 crystals out of the atmosphere that attenuate 
the incident sunlight by Rayleigh scattering. The critical CO2 
partial pressure for the onset of this effect is about 5 to 6 bar. 
However, the cooling e ffect of CO 2 clouds has been challenged 
by Forget & Pierrehum bertI ( Il997l) . CO2 clouds have the addi- 
tional effect of reflecting the outgoing thermal radiation back to 
the surface. The precise inner and outer limits of the climatic 
habitable zone are still unknown due to the limitations of the ex- 
isting climate models. For the present Sun, the HZ is probably 
smaller than between 0.7 to 2 AU, but it is still impossible to give 
a better constraint, particularly for the outer boundary of the HZ. 
For limitations of the planetary habitability of M-type stars see 
iTarteret all (120071) . 

The luminosity and age of the central star play important 
roles in the manifestation of habitability. T he luminosity of Gl 
581 can be obt ained by (1) photometry ( iBonfils et al.l 120051 ; 
lUdrv et al .*20Q7^. and (2) the application of the mass-radius rela- 
tionship (Ribas 2006]) together with the spectroscopi cally deter- 
mined stellar effective temperature of = 3480 K ( Bean et al 
I2OO6I) Both methods yield L = 0.013 + 0.002Ln. IBonfils et al 
(|2005|) consider a stellar age of at least 2 Gyr. 

In the following, we adopt a definition of the HZ previously 
used by franck et al. (2000a b). Here habitability at all times 
does not just depend on the parameters of the central star, but 
also on the properties of the planet. In particular, habitability 
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Fig. 1. Earth system box model. The arrows indicate the differ- 
ent forcing and feedback mechanisms. The bold arrows indicate 
negative feedback operating towards climate stabilization. 



is linked to the photosynthetic activity of the planet, which in 
turn depends on the planetary atmospheric CO2 concentration 
together with the presence of liquid water, and is thus strongly 
influenced by the planetary dynamics. We call this definition the 
photosynthesis-sustaining habitable zone, pHZ. In principle, this 
leads to additional spatial and temporal limitations of habitabil- 
ity, as the pHZ (defined for a specific type of planet) becomes 
narrower with time due to the persistent decrease of the plane- 
tary atmospheric CO2 concentration. 

2. Estimating the habitability of a super-Earth 

2. 1 . Definition of tlie pliotosyntliesis-sustaining iiabitable 
zone 

The climatic habitable zone at a given time for a star with lumi- 
nosity L and effective tem perature differen t from the Sun can 
be calculated according to Jones et a !] (l2006h based on previous 
results by iKasting et al.l { 19931) as 



Lq ■ S out{Te) 



(1) 



with Sin(Te), Sout(Te) givcn as second order polynomials. 

To assess the habitability of a terrestrial planet, an Earth- 
system model is applied to calculate the evolution of the tem- 
perature and atmospheric CO2 concentration. On Earth, the 
carbonate-silicate cycle is the crucial element for a long-term 
homeostasis under increasing solar luminosity. On geological 
time-scales, the deeper parts of the Earth are considerable sinks 
and sources of carbon. 

Our numerical model couples the stellar luminosity L, the 
silicate-rock weathering rate F„r and the global energy balance 
to obtain estimates of the partial pressure of atmospheric carbon 
dioxide Pco^, the mean global surface temperature Tsurf, and the 
biological productivity 11 as a function of time t (Fig. [T]i. The 
main point is the persistent balance between the CO2 sink in the 
atmosphere-ocean system and the metamorphic (plate-tectonic) 
sources. This is expressed through the dimensionless quantities 



where /wr(0 = F^,(t)/F„r,o is the weathering rate, /4(f) = 
Ac{t)/Ac,Q is the continental area, and/si-(f) s S{t)/So is the areal 
spreading rate, which are all normalized by their present values 
of Earth. Eq. (|2| can be rearranged by introducing the geophysi- 
cal forcing ratio, GFR dVolkll 19871) as 



f V/ AT smi,Pc02) - 



fA 



--: GFR(f). 



(3) 



Here we assume that the weathering rate depends only on the 
global surface temperature and the atmospheric CO2 concentra- 
tion. For the investigation of a super-Earth under external forc- 
ing, we adopt a model planet with a prescribed continental area. 
The fraction of continental area relative to the total planetary 
surface /a is varied between 0.1 and 0.9. 

The connection between the stellar parameters and the plane- 
tary climate can be fo rmulated by using a radiation balance equa- 
tion (IWilliamslll998h 



[1 -fl(r,„rf,Pco,)] = 4/«(r,„rf,i^coJ, 



(4) 



where a denotes the planetary albedo, /« the outgoing infrared 
flux, and R the distance from the central star. The Eqs. (O and 
(HJi constitute a set of two coupled equations with two unknowns, 
Tsurf and Pco,, if the parameterisation of the weathering rate, 
the luminosity, the distance to central star and the geophysical 
forcing ratio are specified. Therefore, a numerical solution can 
be attained in a straightforward manner. 

The photosynthesis-sustaining HZ around Gl 581 is defined 
as the spatial domain of all distances R from the central star 
where the biological productivity is greater than zero, i.e.. 



pHZ :={/?! n(Pco,{R, t\ 0) > 0}. 



(5) 



In our model, biological productivity is considered to be solely a 
function of the surface temperature and the CO2 partial pressure 
in the atmosphere. Our parameterisation yields maximum pro- 
ductivity at Tsurt = 50°C and zero pr oductivity for Tsurf < 0°C 
or Tsurf > 100°C or Pco, < 10"^ bar dFranck et al.ll2000ah . The 
inner and outer boundaries of the pHZ do not depend on the de- 
tailed parameterisation of the biological productivity within the 
temperature and pressure tolerance window. Hyperthermophilic 
life forms can tolerate temperatures somewhat above 100°C. 
However, these chemoautotrophic organisms are outside the 
scope of this study. 

2.2. Silicate rock weatliering 

Weathering plays an important role in Earth's climate because 
it provides the main sink for atmospheric carbon dioxide. The 
overall chemical reactions for the weathering process are 



CO2 + CaSi03 
CO2 + MgSi03 



CaCOa + Si02 
MgC03 + Si02 



The total process of weathering embraces (1) the reaction of sil- 
icate minerals with carbon dioxide, (2) the transport of weather- 
ing products, and (3) the deposition of carbonate minerals in sed- 
iments. When combining all these effects, the normalized global 
mean weathering rate /wr can be calculated as 



fwi- — 



flH+,0 



exp 



13.7K 



urf.O 



(6) 



/wr(f) ■ /a(0 = /sr(0, 



(2) 



following IWalker etakl ( Il98ll) . Here the first factor reflects the 
role of the CO2 concentration in the soil, Psoii, with au+ as 
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the activity of in fresh soil-water that depends on Psoii and 
the global mean surface temperature Tsurf . The quantities ah+.o 
and Tsurf.o are the present-day values for the activity and 
the surface temperature, respectively. The activity an* is it- 
self a function of the temperature and the CO2 concentration 
of the soil. The equilibrium constants for the chemical activi- 
ties of the car bon and sulfur systems involved are taken from 
IStumm & Morgan (1981). Note that the sulfur content of the soil 
also contributes to the global weathering rate, but its influence 
does not depend on the temperature. It can be regarded as an 
overall weathering attribute that has to be taken into account for 
the estimation of the present-day value. 

For any given weathering rate, the surface temperature and 
the CO2 concentration of the soil can be calculated in a self- 
consistent manner. Pjoii is assumed to be linearly related to the 
terrestrial biological productivity 11 (see Volk 1987) and the at- 
mospheric CO2 concentration Pco^ . Thus we have 



soil 



^ n / Pco,.o \ ^ Pco, 



(7) 



where Psoii.o, Ho and PcOj.o are again present-day values. 

Since Earth is now harbouring complex life, the weathering 
rates might be lower on a planet without a complex biosphere. 
The parameterisation of biotic enhancement of weathering is 
based only on an increase of the CO2 concentration in the soil. 
This implies a rather weak functional dependence of the weath- 
ering rate on biological productivity. A ten-fold increase in soil 
CO2 concentration relative to the atmosphere only amounts to 
a 1.56 fold increase in the weathering rate due to the present 
biota. This is a significant underestimate, indicating that much 
of the observed biotic amplification of weathering is due to pro- 
cesses other than inc reased soil CO2 concentration. According to 
ISchwartzmanI d 19991) the total amplification due to complex land 
life is at least a factor of 10 and may exceed 100. To explore the 
effect of a strong er biological amplification of weathering a di- 
rect dependence dLenton & von Blohll2001h of weathering rate, 
/wr, on productivity of complex hfe (e.g. land plants), Ilcompiex, 
with amplification factor, auo, can be included: 



1 



ncomplex 1 . 

n ^ |/wr, 



(8) 



where Hcompiex.o is the productivity of the present biosphere. The 
aim of this calculation is to obtain the value for weathering for 
a planet without a complex biosphere. Therefore the weathering 
rate for an Earth with primitive life /„i-,piimitive (Hcompiex = 0) 
is a factor of abio lower than the weathering rate for a complex 
biosphere 



/wr, primitive — /wr i ^bio ^ 1 ■ 

ffbio 



2.3. Thermal evolution model 



(9) 



Parameterised convection models are the simplest models for in- 
vestigating the thermal evolution of terrestrial planets and satel- 
lites. They have been successfully applied to the evolution of 
Mercury, Venus, Earth, Mars, and the Moon dStevenson et al.l 
ll983HSleepl2"00(ih . lFranck & Bounamal(ll995h have investigated 
the thermal and volatile history of Earth and Venus in the frame- 
work of comparative planetology. The internal structure of mas- 
sive terrestrial planets with one to t en Earth masses has been 
investigated bv [Valencia et al.l ( |2006|) to obtain scaling laws for 
total radius, mantle thickness, core size, and average density as 



a function of mass. Similar scaling laws were found for different 
compositions. We will use such scaling laws for mass-dependent 
properties of super -Earths and also mass-inde pendent material 
properties given bv lFranck & Bounamal (Il995h . 

The thermal history and future of a super-Earth has to be de- 
termined to calculate the spreading rate for solving key Eq. (|2]i. 
A parameterised model of whole mantle convection including 
the volatile exchange between the mantle and surface reservoirs 
dFranck & Bounarnall 19951: lFranck|[l998l) is applied. Assuming 
conservation of energy, the average mantle temperature T,„ can 
be obtained as 



^Kpc(Ri-Rl) 



dT, 
dt 



4 

^ = -4nRj„q,„ + -nEiOiRl - rI), (10) 



where p is the density, c is the specific heat at constant pressure, 
qm is the heat flow from the mantle, E(t) is the energy produc- 
tion rate by decay of radiogenic heat sources in the mantle per 
unit volume, and R„, and Rc are the outer and inner radii of the 
mantle, respectively. The radiogenic heat source per unit volume 
is parameterised as 



E(t) = Eoe- 



(11) 



where A is the decay constant of radiogenic heat and the constant 
£0 is obtained from the present heat flux of q,,, - 0.07 Wm"^ for 
an Earth-size planet at 4.6 Gyr. 

The mantle heat flow is parameterised in terms of the 
Rayleigh number Ra as 



Rm - Rr 



Rac 



with 



Ra = 



ga(Tm - TsuYi)(Rm - Rc) 



(12) 



(13) 



where k is the thermal conductivity, Racrit is the critical value of 
Ra for the onset of convection,/? is an empirical constant, g is the 
gravitational acceleration, a is the coefficient of thermal expan- 
sion, K is the thermal diffusivity, and v is the water-dependent 
kinematic viscosity. The viscosity v can be calculated with the 
help of a water fugacity dependent mantle creep rate. It strongly 
depends on the evolution of the mass of mantle water, M,v, and 
the mantle tempera ture, T„j, i.e., v = v{T„„Miy ) and is parame- 
terised according to lFranck & Bounamal dl995h . 

The evolution of the mantle water can be described by a bal- 
ance equation between the regassing flux Freg and outgassing 
flux Fout as 



dM,y 
dt 



iTTiR^n'Ri) 



(14) 



where /bas is the water content in the basalt layer, pbas is the av- 
erage density, t/bas is the average thickness of the basalt layer 
before subduction, S is the areal spreading rate, d„, is the melt- 
ing generation depth and f^- is the outgasssing fraction of water. 
Ryi^o is the regassing ratio of water, i.e., the fraction of subduct- 
ing water that actually enters the deep mantle. The regassing ra- 
tio depends linearly on the mean mantle temperature T„, that is 
derived from the thermal evolution model via 



Rh2o{T,„) = Rt ■ (T,„{Q) - T,n) + Rh.o.^ 



(15) 
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The factor Rj is adjusted to get the correct modern amount of 
surface water (one ocean mass) for an Earth-size planet and 
^HiO.o is fixed at 0.001, i.e., the value is very low at the begin- 
ning of the planetary evolution because of the enhanced loss of 
volatiles resulting from back-arc volcanism at higher tempera- 
tures. 

The areal spreading rate 5 is a function of the average 
mantle temperature T,„, the surface temperature Tsurf, the heat 
flow from the mantle q,„, and the area of ocean basins Aq 
dTurcotte & Schubertlll982h . given as 



q;„nKAa 



(16) 



In order to calculate the spreading rates for a planet with sev- 
eral Earth masses, the planetary parameters have to be adjusted. 
Therefore, we assume 



' M 



0.27 



(17) 



and where Rp is the planetary radius, see [Valencia et al.l (l2006l) . 
The total radius, mantle thickness, core size and average density 
are all functions of mass, with subscript ffi denoting Earth values. 
The exponent of 0.27 has been obtained for super-Earths (M > 
IMe). The values of R,„, Rc, Aq, the density of the planet, and 
the other planetary properties are scaled accordingly. 

The CO2 concentration in the atmosphere Pco, is derived 
from the total mass of atmospheric carbon Catm calculated from 
the balance between sources and sinks according to 



"COj — ~. — *--atm> 

- AnR- He 



(18) 



where yUco, and are the molar weights of CO2 and C, respec- 
tively. The mass dependent pre-factor^//?^, scales with M"*^ *''* w 

M" and has therefore been neglected in our study. 

In Tab. [T]we give a summary of the selected values for the 
parameters used in the thermal evolution model of the 5Me 
and 8M® super-Earth planets. For comparison, the values for 
an Ear th-size planet are also shown. According to [Valencia et al.l 
(l2007h . we assume that a more massive planet is likely to convect 
in a plate tectonic regime similar to Earth. Thus, the more mas- 
sive the planet is, the higher the Rayleigh number that controls 
convection, the thinner the top boundary layer (lithosphere), and 
the faster the convective velocities. In a first order approxima- 
tion, we assume a fixed thickness of the basalt layer and melt- 
ing depth corresponding to relatively low values. Furthermore, 
the initial amount of water M,v(0) scales linearly with the plane- 
tary mass. This might be an underestimate because more massive 
planets tend to accrete more volatiles. 

The geophysical forcing ratios calculated from the thermal 
evolution model for a planet with one, five and eight Earth 
masses for a relative continental area of 0.3 are depicted in 
Fig.|2] It is obvious that the GFR increases with planetary mass 
and is at time zero about 1.8 times higher for SM® and 2.1 times 
higher for 8Me than for the Earth itself. 

2.4. Tidal locking 

According to iPeald (Il977h . the tidal locking radius rj for a 
planet on a circular orbit can be estimated via 



rr= 0.027^:^1 




6 8 10 

Age (Gyr) 

Fig. 2. Geophysical forcing ratios (GFR) for a planet with M = 
IMe (soHd fine), M = SM® (dotted fine) and M - 8Me (dash- 
dotted line) for a relative continental area of 0.3. 



(19) 



where Pq is the original rotation period of the planet, t is the 
time, 2 ' is the dissipation function and Mjtar is the stellar 
mass (all quantities in cgs units). We assume analogously to 
iKastinget al.l (1 1993b Q = 100 and Pq = 13.5 hr 

Planets inside the habitable zone of M stars are tidally 
locked. Due to tidal locking, a weaker intrinsic magnetic field 
is expected. Using simple scaling laws for the planetary mag- 
netic dipole moment m, summarised bV iGrieBmeier et al. (20051). 
we estimate values of about 0.5 m@ and 0.1 me for Gl 581c 
and Gl 58 Id, respectively. The corresponding sizes of ex- 
pected magnetospheres (compressed by coronal winds) can be 
quantifie d by the stando ff distance Rs of the magnetopause 
CKhodachenko et al]|2007l) . The corresponding values for R^ are 
on the order of several planetary radii. Thereby, the surfaces of 
the super-Earth planets Gl 581c and Gl 58 Id are expected to be 
protected from hot coronal winds. 

For these estimates we assume that the super-Earths Gl 581c 
and Gl 581d have at least liquid outer cores. In contrast, accord- 
ing to I Valencia et al. (2006), all super-Earth planets with masses 
higher than 1 M@ have completely solid cores for the preferred 
set of thermodynamic data. Nevertheless, using so-called high 
thermal parameters give warmer interiors and allow super-Earth 
planets to have a completely liquid core. This again points to a 
magnetospheric protection of super-Earth planetary atmospheres 
in Gl 581 from a dense flow of stellar plasma. 

Furthermore, low mass stars have strong XUV irradiations 
during long time periods. Inside the habitable zone, such radia- 
tions can erode the atmosphere of an Earth-size planet and re- 
sult in an addition al limitation on the definition of habitability 
(iGrieBmeier et al.l l2005). Another problem for the habitability of 
tidally locked planets is the freezing out of atmospheric volatiles 
on the dark side of the planet that makes the planets not habit- 
able. Detailed investigations with the help of three-dimensional 
global circulation cUmate models, including the hydrological cy- 
cle, bv jjoshi et al.l ( [19971) and Joshi (2003) showed that approx- 
imately 100 mbars of CO2 are sufficient to prevent atmospheric 
collapse. Therefore, super-Earth planets that generally are as- 
sumed to have more volatiles should not be restricted in their 
habitability by such effects. In case of an eccentric orbit of Gl 
58 Id (e = 0.2), the planet might be trapped in a m : « tidal lock- 
ing (similar to Mercury in the Solar System), thus increasing its 
rotational period. 
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Parameter 




Value 

5Me 


8Me 


Unit 


Description 




5 X 10 ^' 


5 X 10' 


5 X 10' 


m 


average thickness of the basalt layer 


^as 


0.03 


0.03 


0.03 




mass fraction of water in the basalt layer 


/w 


0.194 


0.194 


0.194 




degassing fraction of water 


dm 


40 X 10' 


40 X 10' 


40 X 10' 


m 


melting depth 


k 


4.2 


4.2 


4.2 


T -1 -1 

J s ' m ' K ' 


thermal conductivity 


M„(0) 


4.2 X 10 


2.1 X 10^^ 


5.36 X 10^' 


kg 


initial amount of mantle water 




3,471 X 10^* 


5, 360 X 10' 


6,085 X 10' 


m 


inner radius of the mantle 




6,271 X 10' 


9,684 X 10' 


10, 994 X 10' 


m 


outer radius of the mantle 


T,„m 


3,000 


3,000 


3,000 


K 


initial mantle temperature 


K 


10"^ 


10"* 


10"' 


m- s"' 


thermal diffusivity 


Pbas 


2,950 


4,005 


4,379 


kg m"' 


density of the basalt 


pc 


4.2 X 10 


4.2 X 10 


4.2 X 10 


J m K 


density x specific heat 


Rt 


29.8 X 10"^ 


29.8 X 10"^ 


29.8 X 10"^ 


K-' 


temperature dependence of regassing ratio 


a 


3 X 10-^ 


3 X 10-^ 


3 X 10-5 


K-' 


coefficient of thermal expansion 


p 


0.3 


0.3 


0.3 




empirical constant in Eq. lll2t 


Raciit 


1,100 


1, 100 


1,100 




critical Rayleigh number 


A 


0.34 


0.34 


0.34 


Gyr-' 


decay constant of radiogenic heat 


Eo 


1.46 X 10-^ 


1.46 X 10^^ 


1.46 X 10^^ 


J S-' m-' 


initial heat generation per time and volume 


g 


9.81 


20.6 


25.5 


m s"- 


gravitational acceleration 



3. Results and discussion 

The pHZ around Gl 581 for super-Earths with five and eight 
Earth masses has been calculated for L = 0.011,0.013, and 
O.O15L0. The results for 5M^ are shown in Fig.|3t-c. The simu- 
lations have been carried out for a maximum CO2 pressure of 5 
bar (light colours) and 10 bar (dark colours) neglecting the cool- 
ing effect of CO2 clouds. We assume that the maximum CO2 
pressure of the atmosphere is not limited by the total amount of 
carbon on the planet. Sufficient amounts of carbon are assumed 
to be always available for building up a CO2 atmosphere of up to 
^'max = 10 bar. The biogenic enhancement factor of weathering 
abio has been set to one, i.e., a direct dependence of weather- 
ing on biological productivity according to Eq. ([8]) is neglected. 
Hence, the pHZ is calculated assuming a biotic enhancement of 
weathering similar to the modern biosphere on Earth. The tidal 
locking radius given by Eq. (fT9] ) is also shown. It is evident that 
both planets are well inside the tidal locking radius, assuming 
a stellar age of at least 2 Gyr (iBonfils et al. 1 120051) . The inner 
boundary of the pHZ moves slightly outward, whereas the outer 
boundary decreases nonlinearly with age. Up to a critical age, 
the outer limit is constant and is determined by the maximum 
CO2 atmospheric pressure. Beyond this age, the outer bound- 
ary moves inward due to geodynamic effects. At this point the 
source of carbon released into the atmosphere is too low to pre- 
vent a freezing catastrophe. 

The planet Gl 581c is clearly outside the habitable zone for 
all three luminosities. The luminosity of the central star would 
have to be as low as L = O.OO45L0 to yield habitable solutions 
for this planet. It should be pointed out that Gl 581c is closer 
to its parent star than Venus to the Sun, even with the stellar 
luminosity scaled accordingly. The results for Gl 58 Id are much 
more encouraging (Fig.|4^-c). In particular, for a maximum CO2 
concentration of 10 bar, the planet is inside the pHZ for L > 
0.01 ITLq, a value consistent with the observed luminosity of Gl 
581. Assuming an Earth-like fraction of continental area Uc - 
0.3, this planet would be habitable for a duration of 5.75 Gyr. 
For a water world [Gc = 0.1) the duration is extended to 9 Gyr, 
while for a land world (oc = 0.9) the duration is shortened to 0.3 
Gyr 



Table 2. Maximum outer limit Ro^t of the photosynthesis- 
sustaining habitable zone of Gl 581 for L - O.OOLq as a func- 
tion of the maximum atmospheric CO2 pressure P,nax for three 
different climate models. 
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Due to the tidal locking and the planet's position near the 
outer edge of the pHz, the appearance of complex life is rather 
unlikely on Gl 58Id, i.e., Hcompiex = 0. Life is adversely af- 
fected by low temperatures, small amounts of visible light, and 
insufficient shielding from intense flares. Thus, the simulations 
have been repeated with a parameterisation of weathering ac- 
cording to Eq. ([8]l. The biogenic enhancement factor Cbio has 
been set to 3.6 a ssuming that complex life amplifies weather- 
ing by this factor ( I von Bloh et al.ll2003bh . We find that the dura- 
tion of habitability for primitive Ufe is extended to 8.8 Gyr for 
Qc - 0.3. In gen eral, higher values of gbio extend the life spa n 
of the biosphere ('Lent on & von Blohll200ll:lFranck et al.ll2006l) . 
Incidentallv.iTarter et al.l (l2007h argue that complex life might in 
pri nciple be possible on planets around M stars. 

lUdrv et akl (l2007h determined the orbital eccentricity of Gl 
581c and Gl 58Id as 0.16 + 0.07 and 0.20 + 0.10, respectively. 
Assuming an eccentric orbit does not affect the result that Gl 
581c is outside the pHZ at all times; however, for Gl 58 Id, it 
would imply that it leaves and re-enters the pHZ during its or- 
bital motion for any of the assumed model parameters (i.e., plan- 



D 



13 

3- 




W. von Bloh et al.: The habitability of super-Earths in GUese 581 

(a) 



2 4 6 8 10 12 14 
Age (Gyr) 



1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Relative continental area 

Fig. 3. The pHZ of Gl 581 for a super-Earth (M = 5Me) with a 
relative continental area varied from 0.1 to 0.9 and a fixed stel- 
lar luminosity of (a) 0.01 ILq, (b) O.OOLq, and (c) O.OISLq as 
a function of planetary age. The light colours correspond to a 
maximum CO2 pressure of 5 bar, whereas the dark colours cor- 
respond to 10 bar. For comparison, the positions of Venus, Earth 
and Mars are shown scaled to the luminosity of Gl 58 1 . The light 
grey shaded area denotes the HZ calculated from Eq. ([TJ, while 
the dar k shaded area corresp onds to an extended outer limit fol- 
lowing [Mischnaetal] ( I2OOOI) . The vertical bar at 2 Gyr denotes 
the range of distances due to the (possibly) eccentric orbit. The 
area below the solid black curve is affected by tidal locking. 
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Fig. 4. The pHZ of Gl 581 for a super-Earth (M = 8Me) with a 
relative continental area varied from 0. 1 to 0.9 and a fixed stel- 
lar luminosity of (a) 0.01 ILq, (b) O.OOLq, and (c) O.OI5L0 as 
a function of planetary age. The light colours correspond to a 
maximum CO2 pressure of 5 bar, whereas the dark colours cor- 
respond to 10 bar For comparison, the positions of Venus, Earth 
and Mars are shown scaled to the luminosity of Gl 58 1 . The light 
grey shaded area denotes the HZ calculated from Eq. ([1]), while 
the dar k shaded area corresp onds to an extended outer limit fol- 
lowing |Mischnaeral] ( I2OOOI) . The vertical bar at 2 Gyr denotes 
the range of distances due to the (possibly) eccentric orbit. The 
area below the solid black curve is aff'ected by tidal locking. 



etary climate model and stellar luminosity). But this would not 
thwart planetary habitability since a planet with a sufficiently 
dense atmosphere coul d harbour Ufe even if its orb it is temporar- 
ily outside the HZ, see lWilUams & Pollardl(l2002l) . 



However, it is noteworthy that the possible non-zero eccen- 
tricity values for the Gl 581 super-Earths are highly unlikely 
owing to the sma ll number of radial-velocity measurements by 
IUdrvetalJ(l2007h . As a small number of measurements always 
tends to render an overestimate in the deduced eccentricity, as 
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well as a highly uncertain error bar ( iButleretal J 120061) . a cir- 
cular orb it is strongly preferred in this case, as also done by 
lUdrv et a l. (2007). In this case, the planetary orbit of Gl 58 Id 
is found to stay inside of the pHZ all the time for most of the 
considered planetary climate models and stellar luminosities. 

The ultimate life span of a super-Earth is determined by the 
merging of the inner and outer pHZ boundaries and depends on 
the planetary mass. For a planet older than this ultimate life span, 
no habitability is found. An Earth-like planet with IM® and a rel- 
ative continental area of 0.3 has an ultimate life span of 8.8 Gyr, 
while super-Earth planets with SM® and 8Me have ultimate life 
spans of 11 . 1 Gyr and 11 .9 Gyr, respectively. The critical age and 
ultimate life span is found to decrease with the relative continen- 
tal area It is obvious that an almost completely ocean-covered 
planet ("water world") has the highest likelihood of being hab- 
itable; see also previous models for 47 UMa by iFranck et al.l 
(l2003b . However, for an age of 2 Gyr, habitability for Gl 58 Id 
is not constraint by the outer edge of the pHZ for continental to 
total planetary surface ratios of less than 0.7. In the case of a 
biogenic enhancement factor a^io -3.6 for complex life, habit- 
ability is maintained even for a land world with a continental to 
total planetary surface ratio of 0.9. 

The simulations have been repeated for different maximum 
CO2 pressures Pmax a nd three c limate models (Tab. |2]i. For 
the climate model by IWilUamsl ( Il998l) . Gl 58 Id is habitable 
for f max > 9 bar, for a grey atmosphere model dChamberlainl 
119801) no habit abiUtv can be found, while for the Budyko model 
(lBudvkolfT982h (climate sensitivity 4 K/2 x CO2) habitability is 
attained for f max ^ 4 bar. 

A planet with eight Earth masses has more volatiles than an 
Earth size planet to build up such a dense atmosphere. This pre- 
vents the atmosphere from freezing out due to tidal locking. In 
case of an eccentric orbit of Gl 58 Id (e = 0.2), the planet is hab- 
itable for the entire luminosity range considered in this study, 
even if the maximum CO2 pressure is assumed as low as 5 bar. 
In conclusion, one might expect that life may have originated on 
Gl 58 Id. The appearance of complex life, however, is unlikely 
due to the rather adverse environmental conditions. To get an ul- 
timate answer to the profound question of life on Gl 58 Id, we 
have to await future space missions such as the TPF/Darwin. 
They will a llow for the first tim e to attempt the detection of 
biomarkers jGrenfell et al.ll2007l) in the atmospheres of the two 
super-Earths around Gl 581. 



Kliodachenko, M. L., Ribas, I., Lammer, H., et al. 2007, Astrobiology. 7. 167 
Lenton, T. M., & von Bloh, W. 2001, Geophys. Res. Lett. 28, 1715 
Mischna, M. A., Kasting, J. R, Pavlov, A., & Freedman, R. 2000, Icarus, 145, 
546 

Peale, S. J. 1977, in Planetary Satellites, ed. J. A. Burns, 87 
Ribas, I. 2006. Astrophys. Space Sci.. 304, 89 

Schwartzman, D. W. 1999. Life. Temperature and The Earth: The Self- 
organizing Biosphere. Columbia University Press. New York. 
Sleep, N. H. 2000, J. Geophys. Res., 105, 17563 
Stevenson, D. J., Spohn, T, & Schubert, G. 1983, Icarus, 54, 466 
Stumm, W., & Morgan, J. J. 1981, Aquatic Chemistry, Wiley, New York 
Tarter, J. C, Backus, R R., Mancinelli, R. L., et al. 2007, Astrobiology, 7, 30 
Turcotte, D. L., & Schubert, G. 1982, Geodynamics 
Udry, S., Bonfils, X., Delfosse, X., et al. 2007, A&A, 469, 43 
Valencia, D., O'Connell, R. J., & Sasselov, D. 2006, Icarus, 181, 545 
Valencia, D., Sasselov, D., & O'Connell, R. J. 2007, ApJ, 656, 545 
Volk, T. 1987, Am. J. Sci., 287, 763 

Von Bloh, W., Cuntz, M., Franck, S., & Bounama, C. 2003a, Astrobiology, 3, 
681 

Von Bloh, W., Franck, S., & Bounama, C. 2003b, Geophys. Res. Lett. 30, 1963 
Von Bloh, W., Bounama, C, & Franck, S. 2007, Planet. Space Sci., 55, 651 
Walker, J. C. G., Hays, P B., & Kasting, J. F 1981, J. Geophys. Res., 86, 9776 
Williams, D. M. 1998, The Stability of Habitable Planetary Environments 
Williams, D. M., & Pollard, D. 2002, Int. J. Astrobiology, 1, 61 



References 

Bean, J. L., Benedict, G. R, & Endl, M. 2006, ApJ, 653, L65 

Bonfils, X., Forveille, T, Delfosse, X., et al. 2005, A&A, 443, L15 

Budyko, M. I. 1982, The Earth's Climate: Past and Future 

Buder, R. P, Wright, J. T. Marcy, G. W., et al. 2006, ApJ, 646, 505 

Chamberlain, J. W. 1980, Planet. Space Sci., 28, 1011 

Cuntz, M., von Bloh, W., Bounama, C, & Franck, S. 2003, Icarus, 162, 215 

Forget, F, & Pierrehumbert, R. T. 1997, Science, 278, 1273 

Franck, S., & Bounama, C. 1995, Phys. Earth Planet. Inter, 92, 57 

Franck, S. 1998, Tectonophysics, 291, 9 

Franck, S., Block, A., von Bloh, W., et al. 2000a, Tellus, 52B, 94 
Franck, S., von Bloh, W., Bounama, C, et al. 2000b, J. Geophys. Res., 105 (El), 
1651 

Franck, S., Cuntz, M., von Bloh, W., & Bounama, C. 2003, Int. J. Astrobiology, 
2, 35 

Franck, S., Bounama, C, & von Bloh, W. 2006, Biogeosciences 3, 85 
Grenfell, J. L., GrieBmeier, J.-M., Patzer, B., et al. 2007, Astrobiology, 7, 208 
GrieBmeier, J. M., Stadelmann, A., Motschmann, U., et al. 2005, Astrobiology, 
5,587 

Jones, B. W., Sleep, R N., & Underwood, D. R. 2006, ApJ, 649, 1010 
Joshi, M. M., Haberle, R. M., & Reynolds, R. T. 1997, Icarus, 129, 450 
Joshi, M. M. 2003, Astrobiology, 3, 415 

Kasting, J. F., Whitmire, D. P, & Reynolds, R. T. 1993, Icarus, 101, 108 



